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We study the effects of hole doping on one-dimensional Mott insulators with orbital degrees of 
freedom. We describe the system in terms of a generalized t — J model. At a specific point in 
parameter space the model becomes integrable in analogy to the one-band supersymmetric t — J 
model. We use the Bethe ansatz to derive a set of nonlinear integral equations which allow us to 
study the thermodynamics exactly. Moving away from this special point in parameter space we 
use the density-matrix renormalization group applied to transfer matrices to study the evolution 
of various phases of the undoped system with doping and temperature. Finally, we study a one- 
dimensional version of a realistic model for cubic titanates which includes the anisotropy of the 
orbital sector due to Hund's coupling. We find a transition from a phase with antiferromagnetically 
correlated spins to a phase where the spins are fully ferromagnetically polarized, a strong tendency 
towards phase separation at large Hund's coupling, as well as the possibility of an instability towards 
triplet superconductivity. 



PACS numbers: 71.10.Fd, 05.70.-a, 05.10.Cc 



INTRODUCTION 



In many transition metal oxides, different orbital con- 
figurations are close in energy or even degenerate. Small 
changes in temperature or pressure can therefore lead to 
a complete rearrangement of the electron clouds which in 
turn also strongly influences the magnetic and the trans- 
port properties. Such orbital degrees of freedom play an 
important role, for example, in the manganites, titanates, 
vanadates, and ruthenates. 1,2 Quite common for these 
transition metal oxides is the perovskite crystal structure 
where each transition metal ion is surrounded by an octa- 
hedron of oxygen ions. For a 3d transition metal ion, the 
cubic crystal field then splits the fivefold orbital degen- 
eracy into threefold degenerate t2 9 orbitals and twofold 
degenerate e g orbitals. Here we want to concentrate first 
on the case where we have one electron per site in the t^g 
orbitals with the e g orbitals being inactive (fully occupied 
or empty). Because the onsite Coulomb interactions are 
large the system is a Mott insulator in this case. The 
strongly anisotropic shape of the ti g orbitals means that 
the direction an electron can move to create a virtually 
excited state depends on the orbital it is sitting in. More 
precisely, hopping is only possible between orbitals of the 
same kind, and along a particular crystal axis only two 
out of the three t2 g orbitals are active. This can lead 
to a (dynamical) lowering of the effective dimensionality 
of the system in various ways: Conventional orbital or- 
dering can restrict the hopping to one-dimensional chains 
which can then show typical one-dimensional phenomena 
like a Haldane gap 3 or a Peierls effect. 4 More unconven- 
tional mechanisms like an orbital-driven Peierls effect 5,6 
or spin-orbital nematic states 2 might also render the sys- 
tem quasi one dimensional. In the latter cases, however, 
there will still be a twofold orbital degeneracy. A simple 



Hamiltonian capturing the essential physics is then given 
by 

H = 2JJ2 (SjS j+ x + x) (t 3 t j+1 + y) (1) 



where S is an S — 1/2 spin operator and r a r = 1/2 
orbital pseudospin describing the occupation of the two 
degenerate orbitals active along the chain direction. J = 
4t 2 /U is the magnetic superexchange constant, t the hop- 
ping amplitude, U the onsite Coulomb repulsion, and x, y 
real numbers often treated as free parameters. From a 
microscopic derivation of the effective model (1) it fol- 
lows, however, that x,y are determined by the Hund's 
rule coupling Jh with x = y = 1/4 corresponding to 
Jh = 0. In addition, such a derivation shows that a fi- 
nite Hund's coupling does not only modify x, y but also 
leads to an xxz-type anisotropy of the orbital sector. 2 
This anisotropy is neglected in (1). We will come back 
to the relation between this simple model and more real- 
istic models in section V. 

The spin-orbital model (1) has been intensely 
studied 7-14 and a number of different phases depending 
on x, y have been identified (see e.g. Refs. 11, 15). In 
general, the model has a SU{2) x SU(2) symmetry and 
exhibits an additional Z2 symmetry, interchanging spin 
and orbital degrees of freedom, if x = y. At the special 
point x — y = 1/4 the symmetry is enlarged even further 
to SU(4). This has to do with the fact that at this point 
the Hamiltonian is just a permutation operator of states 
on neighboring sites. The model therefore becomes a ver- 
sion of the Uimin-Sutherland model and is integrable by 
Bethe ansatz. 16 

In this work we want to study the effects of hole doping 
on the spin-orbital model (1). Because states with more 
than one electron per site are effectively forbidden due 
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to the strong Coulomb repulsion U, hole doping of the 
Mott insulator (1) naturally leads us to a generalized t— J 



J 



model 



h.C.^P + 2J^ {(SjSj+l + X) (TjTj+I + 



n 3 n J + l 



(2) 



Here V projects out the doubly occupied states, a is the 
spin, and r the orbital index. As for the one-band t — J 
model it turns out that there is a special point in pa- 
rameter space J/t = 2, x = y = 1/4 where the model 
is integrable by Bethe ansatz. Again, the symmetry is 
enlarged at this point to SU(4\1) (graded SU(5) symme- 
try), the Hamiltonian is a permutation operator of states 
on neighboring sites and falls into the Uimin-Sutherland 
class of models. 

In Sec. II we will investigate the thermodynamics of 
this model at the integrable point with the help of the 
Bethe ansatz and the quantum transfer matrix approach. 
Details of the Bethe ansatz calculation are presented in 
appendix A. In Sec. Ill we briefly introduce the density- 
matrix renormalization group applied to transfer matri- 
ces (TMRG) which we will use to study the thermody- 
namics of model (2) numerically away from the integrable 
point. We will test the accuracy of this method by com- 
paring with exact results at the integrable point. In sec- 
tion IV we will use the TMRG algorithm to study the 
evolution of various phases of the undoped model (1) with 
doping and temperature. In Sec. V we finally consider 
a one-dimensional version of a realistic model for cubic 
titanates which includes the xxz-type anisotropy of the 
orbital sector due to Hund's coupling. We investigate 
the phase transitions, possible tendencies towards phase 
separation as well as superconducting instabilities as a 
function of the strength of Hund's coupling. In Sec. VI 
we present a short summary and our conclusions. 



II. THE INTEGRABLE MODEL 

The integrable SU(4\1) model (2) with J/t = 2, x = 
y = 1/4 has already been studied by Schlottmann. 17 ' 18 
He derived the Bethe ansatz equations and studied the 
ground state properties as well as the elementary excita- 
tions. Kawakami later then derived the critical exponents 
of various correlation functions. 19 Here we want to con- 
centrate on the thermodynamics of this model. Based 
on the quantum transfer matrix approach we derive a set 
of nonlinear integral equations (NLIE) which then are 
evaluated numerically to obtain various thermodynamic 
quantities. Details about the derivation of the NLIE are 
given in appendix A. 

Our results for the thermodynamics in the low tem- 
perature limit can be connected to Schlottmann's and 
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FIG. 1: (color online) Specific heat for the integrable model 
(2) with J/t = 2 and x = y = 1/4 calculated by BA. The 
left inset shows the spin-orbital (v ao ) and the charge velocity 
(v c ) as a function of particle density n. The right inset com- 
pares the BA results (solid lines) with the low temperature 
asymptotics obtained from CFT (dashed lines), see Eq. (4). 



Kawakami's results for the elementary excitations using 
conformal field theory. The SU(4) spin-orbital model, 
i.e., model (2) with n = 1, is known to belong to the 
universality class of the SU(A)i Wess-Zumino-Witten 
(WZW) models, 11,20 so that the central charge c so = 3. 
Similar to the one-band supersymmetric t — J model we 
expect that the critical theory for the hole-doped model 
is a semidirect product of the spin-orbital and the charge 
part so that the free energy at low temperatures is given 
by 



f = eo 



rp2 



(3) 



Here eo is the ground state energy which can be calcu- 
lated by Bethe ansatz 18 , v so (v c ) are the velocities of the 
elementary spin-orbital (charge) excitations, respectively, 
and c c = 1 the central charge of the charge sector. 

In Fig. 1 we show the specific heat as a function of 
temperature for various fillings. According to Eq. (3) the 
specific heat at low temperatures is linear and determined 
by the elementary charge and spin-orbital excitations 
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FIG. 2: (color online) Magnetic susceptibility Xs f° r the inte- 
grable model. The dashed lines denote the zero temperature 
limit according to CFT. 
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FIG. 3: (color online) Compressibility \ c for the integrable 
model. The lines denote the zero temperature limit according 
to CFT. The inset shows the dressed charge as a function of 
density. 



As shown in the left inset of Fig. 1 the velocities of the 
charge and spin-orbital excitations go to zero for n — » 
so that the slope of C diverges in this limit. For n — ► 1, 
on the other hand, only « c -t whereas v so — > Jtt/A. 
This also leads to a diverging slope, however, the charge 
excitations are quickly exhausted so that this behavior 
is only visible at very low temperatures. At higher tem- 
peratures (but still T <C t) the slope then crosses over 
to the value given by the spin-orbital excitations only. 
Finally, for n = 1, we have C = 2T in the whole confor- 
mal regime. In the low temperature limit, the BA results 
which we obtained by the quantum transfer matrix ap- 
proach indeed agree perfectly with the CFT result (4) 
using the velocities determined according to Ref. 18 (see 
right inset of Fig. 1). 

The magnetic susceptibility Xs as a function of tem- 
perature is shown in Fig. 2. From CFT we expect 
Xs = 1/iTVso at zero temperature. Because v so vanishes 
for n — + 0, the magnetic susceptibility diverges in this 
limit. For n = 1, on the other hand, we have v so = Jit/A 
so that Xs — 2/7T 2 . Note, however, that logarithmic 
corrections are expected at low temperatures similar to 
the Heisenberg chain. Therefore the susceptibility will 
approach the zero temperature limit predicted by CFT 
with infinite slope. This explains why even at the lowest 
temperatures shown in Fig. 2 the susceptibility curves 
obtained by BA still deviate significantly from the zero 
temperature limit. In Fig. 3 we show the compressibility 
Xc for various densities. For the compressibility CFT pre- 
dicts Xc = K c /(ttv c ) where K c is the Luttinger parameter 
of the charge sector. Using the BA we can calculate the 
so called dressed charge 19 ' 21,22 , £ C (Q), which is related to 
the Luttinger parameter by K c — £ C (Q) 2 - The dressed 
charge as a function of density is shown in the inset of 
Fig. 3. Because the charge velocity vanishes for n — > 
and n — > 1 (see Fig. 1) we have a diverging compress- 
ibility in both limits. In addition, we notice that even 
at intermediate densities the compressibility is large in- 
dicating that the integrable point is not that far from 




FIG. 4: (color online) Drude weight D for the supersymmetric 
point J = 2t (black solid line) and for J = (red dashed 
line), respectively. In the insets the regions of small and large 
densities are shown in detail. 



a phase separated state. In general, phase separation is 
promoted by an increasing ratio of J/t as is well known 
in the SU (2) t — J model. 23,24 We will come back to this 
point in Sec. IV. 

Finally, we want to consider the transport properties of 
the integrable model at zero temperature. We can write 
the real part of the conductivity at zero momentum as 
a'(uj) — 2ttDS(lo) +<7 r eg(w) where D is the Drude weight 
and crreg^) the regular part. With the help of conformal 
field theory and Bcthe ansatz we find 



D = 



K c v c 
~2^~ 



2tt 



(5) 



Thus, the Drude weight follows directly from the dressed 
charge shown in the inset Fig. 3 and the charge velocity 
shown in the left inset of Fig. 1 and is depicted in Fig. 4. 
D vanishes for n — > where the number of charge carriers 
vanishes, and for n — > 1 where the Mott gap opens up. 
Quite surprisingly, the Drude weight for the supersym- 
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metric point appears to be symmetric around n = 1/2 
although the Hamiltonian (2) does not possess such a 
symmetry. To understand this behavior it is instructive 
to study in addition the Drude weight for J = 0. In this 
case only the kinetic energy part of (2) remains. Then 
the spin and orbital indices do not matter and due to the 
projection operators the model becomes equivalent to a 
one-band, spinless fermion model. Note that the Mott 
gap at n = 1 is therefore still incorporated. The Drude 
weight can now be calculated from 



e(*) - e(0) = D&/N 2 + 0(N- 3 ) 



(6) 



with e being the ground state energy per site, and $ 
a field describing a twist in the boundary conditions, 

Cj+N 



We find 



t 



D(J = 0) = — sin 7m 

7T 



(7) 



which is shown as dashed curve in Fig. 4. Now the Drude 
weight is indeed symmetric around n = 1/2 because the 
spinless fermion model is particle-hole symmetric. The 
comparison with the supcrsymmctric case near n ~ in 
the left inset of Fig. 4 shows that the Drude weights coin- 
cide in this limit. This is expected because the exchange 
interaction J becomes irrelevant in the dilute limit. For 
n ~ 1 shown in the right inset of Fig. 4, on the other 
hand, we see that the two curves do not coincide. In this 
limit J cannot be neglected. This means that D for the 
supersymmetric case is not symmetric around n = 1/2 
but the deviations from this symmetry are very small. 
Similar observations have been made for the one-band 
supersymmetric t — J model by Kawakami and Yang, 
Ref. 25. 



III. DENSITY-MATRIX RENORMALIZATION 
GROUP 

The density-matrix renormalization group applied to 
transfer matrices (TMRG) is based on a mapping of a 
one-dimensional quantum system to a two-dimensional 
classical one by means of a Trotter-Suzuki decomposition. 
In the classical model one direction is spatial whereas 
the other corresponds to the inverse temperature. For 
the classical system a so called quantum transfer matrix 
(QTM) is defined which evolves along the spatial direc- 
tion. At any non-zero temperature the QTM has the cru- 
cial property that its largest eigenvalue A is separated 
from the other eigenvalues by a finite gap. The parti- 
tion function of the system in the thermodynamic limit 
is therefore determined by Ao only, allowing it to perform 
this limit exactly. The Trotter-Suzuki decomposition is 
discrete so that the transfer matrix has a finite number 
of sites or local Boltzmann weights M. The temperature 
is given by T ~ (eM) _1 where e is the discretization pa- 
rameter used in the Trotter-Suzuki decomposition. The 
algorithm starts at some high-temperature value where 



M is so small that the QTM can be diagonalized exactly. 
Using a standard infinite-size DMRG algorithm, sites are 
then added to the QTM leading to a successive lowering 
of the temperature. The TMRG algorithm is described in 
detail in Refs. 26, 27, 28 and has been applied to a num- 
ber of one-dimensional systems such as frustrated and 
dimerized spin chains, 29 the Kondo lattice model, 30 the 
t — J chain 28 ' 31 and ladder 32 as well as to the extended 
Hubbard model. 33 

To obtain insight into the physical properties of the 
doped spin-orbital model we will, in particular, be inter- 
ested in the behavior of two-point correlation functions. 
At finite temperatures we expect that any two-point cor- 
relation function of a local operator 0(r) decays expo- 
nentially with distance r 



(O(l)O(r)) - (0(l)){0(r)) =J2M n e 



(8) 



Here M n is a matrix element, the correlation length, 
and k n the corresponding wave vector. Note, that in 
the asymptotic expansion (8) infinitely many correlation 
lengths appear. Within the TMRG algorithm, correla- 
tion lengths and corresponding wave vectors are deter- 
mined by next-leading eigenvalues A„ of the QTM 



C 1 = In 



Ao 
A„ 



k n = arg 



Ao 



(9) 



The long-distance behavior of the correlation function is 
then dominated by the correlation length £ Q belonging 
to the largest eigenvalue A a (a ^ 0) with M a ^ 0. 

Apart from spin-spin (orbital-orbital) (S(l)S(r)) 
((r(l)r(r))) two-point correlation functions we are also 
interested in pair correlation functions to investigate pos- 
sible superconducting instabilities. For the model (2) we 
can define the following singlet and triplet pair correla- 
tion functions 



G tt (r) 
G ss (r) 
G ts (r) 
G st (r) 



(c Ta (r + l)c TQ (r) C } a (2) C t Ta (l)) 
<c To (r + l) Ci6 (r)c| (2) C { 6 (l)> 
(c Ta (r + l) CTb (r)c{ a (2) C } b (l)) 
(c Ta (r + l)c ia (r) C } a (2)4 Q (l)) . 



(10) 



Here a, b denote the r z -componcnt and f, I the S z - 
componcnt. For each of these pair correlations an asymp- 
totic expansion (8) exists and relation (9) can be used to 
numerically determine the corresponding leading corre- 
lation lengths. 

The discrete Trotter parameter e leads to a systematic 
error in the free energy of order e 2 . In the calculations 
presented here we have chosen e = 0.05 so that this error 
is expected to be of the order 10~ 3 — 10~ 4 . More impor- 
tant is the error due to the truncation of the Hilbert space 
in each DMRG step. This error is difficult to estimate 
but accumulates with each DMRG step finally leading to 
a breakdown of the numerics at low temperatures. In the 
calculation presented here we will keep N = 240 — 360 
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FIG. 5: (color online) (a): Free energy / calculated by TMRG 
(symbols) with N = 360 states kept compared to the exact 
solution (lines). The black circles (black solid lines) denote 
the result for n — —1.7, the red squares (red dashed lines) 
for fi — —1.9, and the blue diamonds (blue dot-dashed lines) 
for fi = —2.5, respectively, (b): Absolute error |A/| of the 
TMRG results presented in (a), (c): TMRG results for the 
density compared to the exact solution with symbols and lines 
denoting the same chemical potentials as in (a), (d): Absolute 
error |A n | of the TMRG results in (c). 



states as basis for the truncated Hilbert space. To decide 
down to which temperatures the TMRG is reliable we 
show results for the free energy and density at the inte- 
grable point in Fig. 5. It is important to note, that we 
perform the numerical calculations in a grand canonical 
ensemble, i.e., we fix the chemical potential and not the 
particle density. In particular for small doping levels it 
is, however, possible to find a chemical potential so that 
the density depends only very weakly on temperature as 
shown in Fig. 5(c). Note that for T — > oo we always have 
n — > 4/5 because we have 5 states locally with one state 
corresponding to the empty site. The absolute errors in 
the free energies and in the densities stay smaller than 
10~ 2 for temperatures down to T/t ~ 0.1 as shown in 
Fig. 5(b) and (d), respectively. This accuracy is com- 
pletely sufficient to study the thermodynamic properties 
of model (2), and temperatures of the order T/t ~ 0.1 
are low enough to identify the ground state as well. 



IV. THE ij-MODEL 

In this section we want to investigate model (2) away 
from the integrable point. First, we want to demonstrate 
that the integrable point is indeed already close to a state 
with phase separation. While keeping x = y = 1/4 we 
now set J/t — 3 and show in Fig. 6 the density as a func- 
tion of chemical potential for various temperatures. For 
T/t — > the density is zero for /i < —2.73 and equal to 
one otherwise. This means that the compressibility cal- 



FIG. 6: (color online) The density n as a function of chemical 
potential fi for different temperatures. Here J/t = 3 and 
x = y — 1/4. The lines are guides to the eye. 



culated at any fixed density is divergent and the ground 
state therefore phase separated. The reason for phase 
separation is obvious: For large J/t the system tries to 
maximize its magnetic exchange energy which is achieved 
by separating the particles from the holes. 

For all transition metal oxides we expect in general 
J = At 2 /U < t which is equivalent to At < U. Exact 
values for J/t depend on the considered compound but 
values of J/t ~ 0.3 — 0.5 are typical. As a representative 
value we will concentrate in the following on J/t = 0.5. 



A. The dimerized phase 

For x = y = 1/2 the undoped model is in a dimerized 
phase and spin and orbital excitations are gapped. 11,14 
Spin-Peierls-type instabilities are a generic feature of sys- 
tems with coupled spin and orbital degrees of freedom 
and have been investigated in more detail in Refs. 6, 
34. Here we want to study how dimer order and exci- 
tation gap evolve with doping. In Fig. 7 the magnetic 
susceptibility for various doping levels is shown. Note, 
that due to the Zi symmetry spin and orbital sectors are 
equivalent and spin and orbital susceptibility therefore 
identical. In the undoped case a spin gap A is clearly 
visible. In Ref. 14 this gap has been found to be of the 
order A = 0.090 ±0.005. With increasing hole concentra- 
tion the spin gap becomes smaller but is still detectable 
numerically at a chemical potential (i = 1.0 which corre- 
sponds to a density n = 0.85 at low temperatures. For 
larger doping levels the spin excitations seem to become 
gapless so that the system apparently turns into a Lut- 
tinger liquid. The long-range dimer order, on the other 
hand, seems to break down immediately when holes are 
added to the system. In the case of algebraically de- 
caying correlations at zero temperature we expect the 
corresponding correlation length to diverge as £ ~ 1/T 
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FIG. 7: (color online) Spin (orbital) susceptibility for J/t — 
0.5 and x = y = 1/2 and different chemical potentials. Insets: 
(a) Corresponding densities as function of temperature, (b) 
Leading spin (orbital) dimer correlation length £d. 



whereas £ will diverge stronger than 1/T in the case of 
true long-range order. In inset (b) of Fig. 7 we therefore 
show the leading spin (orbital) dimer correlation length 
£d multiplied by T. In the undoped case the diver- 
gence of T£d indicates that the ground state has indeed 
long-range dimer order, whereas for a chemical potential 
/i = 1.1, corresponding to n ~ 0.9 at low temperatures, 
T£d is decreasing with temperature. Here we expect £u 
to stay finite so that T£d — > for T — » in accord with 
the numerical data. 



B . Ferromagnetic / Antiferromagnetic-phase 

Here we want to consider J/t — 0.5 with x = 0.5, 
y = —0.5 where the undoped model shows ferromag- 
netism in the spin sector and algebraically decaying an- 
tiferromagnetic correlations in the orbital sector. 11 As 
shown in inset (a) of Fig. 8 we find that the density for a 
chemical potential fi — —0.1 is almost constant, n ~ 0.84, 
for temperatures T s [0,2]. This allows us to compare 
the undoped model directly with this slightly doped case. 
In the main figure the spin susceptibility Xs as a function 
of temperature is shown. As expected, \s becomes sup- 
pressed with doping but still diverges for T — > indicat- 
ing long-range ferromagnetic order in both cases. In inset 
(b) of Fig. 8 the nearest-neighbor spin-spin and orbital- 
orbital expectation values are shown. In the undoped 
case (SiSi+x) — > 1/4 for T — > as expected for ferro- 
magnetic order. When the spins order ferromagnetically, 
then, according to Hamiltonian (1), we have an effec- 
tive antifcrromagnetic coupling for the orbitals so that 
(TiT i+1 ) -» -In. 2 + 1/4 for T -» 0. This is the value 
for an antiferromagnetic Heisenberg chain known from 
Bethe ansatz. In the slightly doped case both the ferro- 
magnetic spin and the antiferromagnetic orbital correla- 



FIG. 8: (color online) TMRG results for J/t = 0.5 with x = 
1/2 and y = —1/2. Main figure: Spin susceptibility in the 
undoped case and for fj, — —0.1 as function of temperature. 
Insets: (a) For y, — —0.1 the density n ~ 0.84 is almost 
independent of temperature, (b) Nearest-neighbor spin-spin 
and orbital-orbital expectation values for the undoped model 
and for ju = -0.1 (n ~ 0.84). The circles on the T = 
axis denote 1/4 and — In 2 + 1/4, respectively, (c) Charge 
compressibility for y, — —0.1. 
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FIG. 9: (color online) Leading correlation lengths for J/t = 
0.5, x — 1/2, y — —1/2. The leading spin-spin correlation 
length is non-oscillating (k = 0) whereas the orbital-orbital 
and the pair correlations show incommensurate oscillations at 
low temperatures as depicted in the inset. 



tions become weaker. The charge compressibility shown 
in inset (c) of Fig. 8 is nonzero for T — > indicating that 
the charge excitations are gapless. 

This leads us to the question if the system has alge- 
braically decaying pair correlations and if so, which one 
of the pair correlations defined in (10) dominates. In 
Fig. 9 we show some of the leading correlation lengths. 
Dominant is the spin-spin correlation length with wave 
vector k = which diverges as 1 /T indicative of the fer- 
romagnetic order in the ground state. 35 Next, we find a 
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correlation length which belongs to the asymptotic ex- 
pansion of the orbital-orbital correlation function. Inter- 
estingly, the associated wave vector is given by k = n for 
T/t > 0.28 but becomes incommensurate at lower tem- 
peratures. At T/t ~ 0.28 the correlation length shows 
a cusp. The explanation for this cusp is as follows: In 
the asymptotic expansion (8) of the orbital-orbital cor- 
relation function we have two correlation lengths with 
wave vector k — it. At T/t ~ 0.28 these two correla- 
tion lengths cross and the oscillations of the now leading 
correlation length become incommensurate. For T — > 
we find that k -> ±2k F = ±irn w ±2.64 for n ~ 0.84. 
So at T = the incommensurate oscillations just re- 
flect the incommensurate filling of the system. Similar 
crossover phenomena in the leading correlation length at 
finite temperature have also been observed in the xxz 
model in a magnetic field 36 and in the one-band t — J 
model at incommensurate filling. 28 ' 31 

For the pair correlations defined in (10) we find that 
the spin triplet/orbital singlet correlation Gts has the 
largest correlation length £ ts . At low temperatures we 
find £ts ~ indicating that G ts will decay alge- 

braically at zero temperature. The associated wave vec- 
tor is given by kt. s = for T/t > 0.5 but becomes incom- 
mensurate at lower temperatures. The numerical data 
seem to indicate that for T — > the oscillations become 
commensurate again with k ts — w. In the ground state 
this would mean that G ts (r) ~ (— l) r /r x with some crit- 
ical exponent x. 

In the ferromagnetic/antiferromagnetic-phase consid- 
ered in this section a coupling between spin-orbital chains 
might therefore induce true long-range triplet supercon- 
ductivity. Here triplet superconductivity arises from the 
coupling of the spins with the orbital pseudospins. The 
degenerate orbitals order antiferromagnetically leading 
to an effective ferromagnetic coupling for the spins. 

V. A ONE-DIMENSIONAL VERSION OF A 
REALISTIC MODEL FOR CUBIC TITANATES 

In recent years a lot of interest has focused on the Mott 
insulator LaTiC>3. 2 ' 37 Here the octahedron of oxygen ions 



surrounding each Ti 3+ is nearly perfect. This opens up 
the possibility that the orbital degeneracy is not lifted by 
lattice distortions and that the orbitals act as additional 
quantum degrees of freedom. Starting from the ideal case 
of completely degenerate i 29 -orbitals one can derive a su- 
perexchange model similar to (1). Here the two orbitals 
represented by the orbital pseudospin r depend on the 
bond direction. Along the c-axis, for example, only the 
tig levels of xz and yz symmetry are active and repre- 
sented by r, whereas t stands for the xy and xz orbital 
if the bond is along the a-axis. Having different orbital 
pairs active along each spatial direction necessarily frus- 
trates the one-dimensional physics discussed in the previ- 
ous sections of this paper and might lead to a liquid state 
with short range Sl/(4)-type correlations. 37 Nevertheless, 
a directional, nematic state where the system makes full 
use of the orbital quantum fluctuations say along the c- 
axis with active xz and yz orbitals while the xy orbital is 
empty, thus preventing orbital fluctuations in the other 
two directions, might be close in energy and could possi- 
bly be realized in LaTi0 3 under pressure. 2 An example, 
where such orbital selection leading to a strongly direc- 
tional spin-orbital state probably happens is YVO3. In 
this system, however, we have an effective spin S = l. 5,6 



It is important to take the Hund's rule splitting of 
the virtual excited states into account when deriving the 
superexchange Hamiltonian for LaTiC>3. This leads, in 
particular, to an xxz-type anisotropy of the orbital sec- 
tor, i.e., contrary to (1) the Hamiltonian no longer has a 
SU(2) x SU(2) symmetry. If we consider the hole-doped 
case and add the hopping of the holes to the superex- 
change Hamiltonian for LaTiC>3 given in Ref. 2, we obtain 
again a t — J type Hamiltonian which can be represented 
as 



3 a , T 

+ 2J eif [{SjS j+1 + x) (tj-Tj+i + 5t/t/ +1 + yn jnj+1 ) - -n n 3+x + tt/t/ +1 



(11) 



with parameters If we ignore the splitting of the virtually excited states 

due to Hund's coupling we have r\ = r 2 = r% = 1 so that 

j _ J_i r , \ x _ 1 , 1 ri ~ r2 the Hamiltonian (12) is equivalent to the Hamiltonian 

e ° 2 ' 4 2 n + r 2 (2) with x = y = 1/4 and J eff = J. For finite Hund's 

1 1 n - r 2 1 r 3 - r 2 r 2 r 3 - r 2 

y = 1 , = 

4 2 n + r 2 6 n + r 2 3 n + r 2 

2 n(5r 2 + r 3 ) 2 n(r 2 - r 3 ) 

3 (rx+r 2 ) 2 ' 7 3 (n+r,) 2 ' 1 ' 



8 




ol ' i ' lJ ' ±= 

0.1 ^ 0.2 0.3 



(b) 


X__________ 1 ~~~ 1 




1 1 

y 1 — — ■ 




1, 7~~~~~~~~Tt^ 



I i I i U i u_z 

0.1 0.2 0.3 



FIG. 10: (color online) (a) Parameters J e g for J = 0.5 and 
z as denned in Eq. (12) as a function of r\ = Jh/U. (b) 
Parameters x and y as a function of rj. (c) Parameters 8 
and 7 related to the icicz-anisotropy of the orbital sector as a 
function of rj. The dashed blue lines in (a) and (b) indicate 
the approximate values for the phase transitions described in 
the text. 

coupling Jh we have 

1 1 1 , . 

n = 5- , r 2 = , r 3 = 13) 

where 77 — Jh/U. 

From optical data and first principle calculations for 
LaTi03 one finds the approximate values for the hop- 
ping amplitude t ~ 0.3 eV, the onsite Coulomb repulsion 
U ~ 2.8 eV, and the Hund's rule coupling J H ~ 0.6 eV. 2 
For the model (11) this means that J = 4t 2 /U ~ 0.13 
eV, J/t = At/U ~ 0.43, and rj = J H /U ~ 0.21. It is 
therefore reasonable to set again J/t = 0.5 as in the pre- 
vious section. The parameters (12) as a function of 77 for 
this value of J/t are shown in Fig. 10. Although 6 and 
7 become nonzero for finite 77 thus destroying the SU (2) 
symmetry of the orbital sector, their values remain small 
in the physical regime for 77 depicted in Fig. 10. It is 
therefore indeed reasonable to neglect this anisotropy in 
a first approximation as has been done in the previous 
section. Furthermore, we also find that z ~ 1 at least 
up to 77 ~ 0.2 so that the variation in z can also be 
neglected. We are then back to Hamiltonian (2) with 
J replaced by an effective superexchange scale J e g and 
with x and y being functions of the single parameter 77 
only. From this observation we can infer the basic prop- 
erties of this model: For 77 < 0.2 the model will be in 
a "reseated SU(4) phase", i.e., a phase where the same 
field theory as at the SU(4) symmetric point describes 
the low-energy properties but with spin and orbital veloc- 
ities which are rescaled and no longer equivalent. Strictly 
speaking this is only correct without orbital anisotropy 
(5 = 7 = 0). Depending on the spin order, the orbital 
anisotropy might become Ising-like so that the orbital 




FIG. 11: (color online) Density as a function of the chemical 
potential for J/t = 0.5 and rj = 0.3. The inset shows the same 
for 77 = 0.25. The lines are guides to the eye. 



excitations become gapped. However, even if this hap- 
pens the orbital gap caused by this mechanism will be 
extremely small. For 77 > 0.2 we expect to enter a phase 
with ferromagnetically ordered spins and antiferromag- 
netic correlations in the orbital sector. At the same time 
J e ff increases with increasing 77 so that we expect a phase 
separated state if 77 > 0.3 in the doped case. 

In the following, we present numerical results for the 
full model (11) with J/t = 0.5 and the parameters as 
given in Eqs. (12,13). In Fig. 11 the density as a function 
of chemical potential is shown for 77 = 0.3 and 77 = 0.25. 
As in Fig. 6 we see that the density for 77 = 0.3 and 
T/t — > jumps from zero to one. Here the jump occurs 
at a chemical potential fi w —3.5. Again this indicates a 
diverging compressibility in a canonical ensemble for all 
densities and confirms the expected phase separation at 
large 77. For 77 = 0.25 (shown the in the inset of Fig. 11), 
on the other hand, the ground state is not phase sepa- 
rated. 

Next, we consider the slightly doped case, n ~ 0.8 — 
0.85, for different parameters 77. The nearest-neighbor 
correlation functions (SjSj+i) and (TjTj+i) presented 
in Fig. 12 show clearly that the spin correlations are an- 
tiferromagnetic for 77 = 0.1 and ferromagnetic for 77 = 0.2 
and 77 = 0.25. The orbital correlations, on the other 
hand, are antiferromagnetic in all three cases. To fix 
the critical value for 77 where the phase transition occurs, 
we consider in Fig. 13 the nearest-neighbor spin corre- 
lation for the undoped model (n = 1) as a function of 
temperature T and Hund's coupling 77. The data show 
that the spins in the ground state are fully polarized if 
77 > ?7 C w 0.2. For r\ < r\ c the spin correlations are an- 
tiferromagnetic in the ground state. In this case, how- 
ever, (SjSj+i) can be non-monotonic as a function of 
temperature and even larger than zero in a certain tem- 
perature range if 77 is close to rj c . The expectation value 
for {SjSj+i) jumps at zero temperature from some 77- 
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FIG. 12: (color online) TMRG results for model (11) with 
J = 0.5 and r\ — 0.1, 0.2,0.25. The chemical potential is set 
to fi = 0.6 for r) = 0.1, fi = 0.0 for ij = 0.2, and /i = -0.8 for 
r\ = 0.25, respectively, (a) Nearest-neighbor spin-spin corre- 
lation {SjSj+i), (b) nearest-neighbor orbital-orbital correla- 
tion (TjTj+i), and (c) the density as a function of tempera- 
ture. Note, that for the chemical potentials chosen here the 
densities depend only weakly on temperature, n ~ 0.8 — 0.85. 
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FIG. 14: (color online) Leading correlation lengths for J = 
0.5, r) = 0.25 and a chemical potential n — —0.8 correspond- 
ing to a density n ~ 0.84 at low temperatures (see Fig. 12(c)). 
The orbital-orbital and the density-density correlations show 
incommensurate oscillations at low temperatures as depicted 
in the inset whereas the spin-spin and the spin triplet /orbital 
singlet correlations are both non-oscillating. 



v 




FIG. 13: (color online) TMRG results for the nearest- 
neighbor correlation function (SjSj+i) in the undoped case 
for J = 0.5 and different rj. The red dashed curve denotes the 
result for r] — 0.2. The spins are fully polarized in the ground 
state, (SjSj+i) = 1/4 (black dot), if r? > 0.2. 



dependent value for rj < rj c to 1/4 for r/ > rj c . The phase 
transition driven by r\ is therefore first order. Very inter- 
estingly, the phase transition occurs at a value for rj which 
is close to the one expected for LaTiC>3. In addition to 
the frustration of the orbital sector in this compound due 
to different orbital pairs being active along each direction, 
this closeness to the phase transition might also be im- 
portant to understand the peculiar physics of LaTiC>3 . It 
might, in particular, be a contributing factor to the small- 
ness of the ordered moment in the G-type antiferromag- 
netic structure. 38 Note, however, that also other factors 
like the Ti-O-Ti bond angle are very important in deter- 



mining whether the spins order ferro- or antiferromag- 
netically. The different magnetic properties of YTiO\3 
(ferromagnetic spin order) and LaTiC>3 (G-type antifer- 
romagnetic order), for example, have been ascribed to a 
small variation in this angle. 39 Nevertheless, this varia- 
tion in bond angle can be mimicked to a certain degree 
by increasing rj so that the phase transition situated at 
t] c rs 0.2 in the one-dimensional model is indeed impor- 
tant to understand the physics of the cubic titanates. 

Finally, we again want to study possible pairing insta- 
bilities in the slightly doped case. Here we concentrate 
on rj = 0.25 with J = 0.5 and fi — —0.8 as in Fig. 12. 
Results for the leading correlation lengths in this case are 
presented in Fig. 14. We find that the spin-spin, orbital- 
orbital, density-density, as well as the spin triplet/orbital 
singlet correlation lengths all diverge as l/T for T — > 
indicating that these correlations will decay algebraically 
at zero temperature. As in Sec. IV B we find that the os- 
cillations of the orbital-orbital and the density-density 
correlation become incommensurate in the low temper- 
ature regime. Again k — > ±2kp — ±irn ±2.64 for 
n ~ 0.84 reflecting the incommensurate filling of the sys- 
tem. The spin triplet/orbital singlet correlation length 
is larger than any of the other correlation lengths as- 
sociated with the pair correlations defined in (10). If 
superconductivity can be stabilized at all in a possible 
nematic phase of LaTi03 or YTi03 there is therefore the 
possibility that it will be of triplet character. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we studied the thermodynamic proper- 
ties of hole-doped one-dimensional Mott insulators with 
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orbital degrees of freedom. We described such systems 
in terms of generalized (multi-band) t — J models. Ne- 
glecting the Hund's rule splitting of the virtually excited 
states we were led to a model which is integrablc for one 
specific value of the ratio J/t. At this point the model 
becomes SU(4\1) symmetric (graded SU(5) symmetry) 
and belongs to the so called Uimin-Sutherland class of 
models. 16 The integrability at this particular point is 
analogous to the integrability of the usual t — J model 
at the supersymmetric point. Ground-state properties 
of the 5C/(4|1) symmetric spin-orbital model have been 
first investigated using Bcthe ansatz by Schlottmann and 
Kawakami. 17-19 Here we presented a set of nonlinear in- 
tegral equations also based on the Bethe ansatz which 
allowed us to study the thermodynamics. Using confor- 
mal field theory we have been able to connect our new 
results for the thermodynamics at low temperatures with 
the results for the ground state and the elementary exci- 
tations obtained by Schlottmann and Kawakami. 

In the second part of the paper we used the density ma- 
trix renormalization group applied to transfer matrices 
(TMRG) to study the thermodynamics of the two-band 
t — J model away from the integrable point. By compar- 
ing with Bcthe ansatz results at the integrable point we 
first demonstrated that the obtained numerical results 
are accurate down to low temperatures T/t ~ 0.05. For 
large values of J/t we then showed that the ground state 
becomes phase separated. Next, we studied the effects of 
hole doping on various phases of the undoped model. If 
the ground state of the undoped model is dimerized then 
the associated spin gap persists up to relatively large hole 
concentrations. The long-range nature of the dimer or- 
der, however, seems to break down immediately upon 
hole doping. Starting from the phase of the undoped 
model where the spins are fully ferromagnetically polar- 
ized and the orbitals show antifcrromagnetic correlations 
we found that upon hole doping the spin triplet/orbital 
singlet pair correlation dominates among the various pos- 
sible pair correlation functions. This correlation function 
will decay algebraically at zero temperature so that in- 
terchain couplings might stabilize true triplet supercon- 
ductivity in this phase. 

In the last part we used the TMRG algorithm to study 
the effects of hole doping for a one-dimensional version of 
a realistic model for cubic titanates. The model can be 
written in a form making it very similar to the one con- 
sidered earlier. Previously independent parameters, how- 
ever, now become functions of the ratio 77 of Hund's cou- 
pling Jh and the onsite Coulomb repulsion U, and there 
is an additional xxz-type anisotropy of the orbital sec- 
tor. The effective superexchange coupling is a function of 
77 as well and increases monotonically with increasing 77. 
For 77 > 0.3 this leads to a phase separated ground state. 
In addition, we find a phase transition at ?7 C w 0.2 be- 
tween a state with antiferromagnetically correlated spins 
(77 < T] c ) and a state with fully ferromagnetically polar- 
ized spins (77 > r] c ). The phase transition is first order. 
For 77 > rj c the dominating pair correlation is again of 



spin triplet/orbital singlet character and might lead to 
a true superconducting instability if interchain couplings 
are present. Interestingly, realistic 77-values for LaTi03 
are close to r\ c . Although a purely one-dimensional model 
is not appropriate for this compound, the nature of lo- 
cal correlations might be correctly captured and 77 being 
close to ?7 C might contribute to the peculiar properties 
of this compound, in particular, to the extremely small 
G-type magnetic moment. 
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APPENDIX A: NONLINEAR INTEGRAL 
EQUATIONS 

The integrable SU(4\1) model admits the calcula- 
tion of exact results for the thermodynamics. Here the 
largest eigenvalue of the QTM can be obtained by Bethe 
ansatz. 40 The number of Bcthe ansatz equations, how- 
ever, diverges in the limit M — > 00. It is thus necessary 
to encode the Bethe ansatz equations into an alternative 
form for which the limit can be taken analytically. This 
can be done by defining suitable auxiliary functions in the 
spirit of Refs. 41,42, which are shown to be determined 
by a closed set of only finitely many coupled nonlinear 
integral equations (NLIEs). 

The rigorous derivation depends on the explicit knowl- 
edge of the auxiliary functions in terms of the Bcthe 
ansatz roots for finite M. Unfortunately these are un- 
known for the SU(4\1) model. Yet we are able to con- 
jecture the complete set of coupled NLIEs by general- 
izing the structure that has been found for two closely 
related models in the Uimin- Sutherland class, namely 
the SU(2\1) and the SU(4) models. 43 - 44 For the SU(4\1) 
model we thus expect a total number of 15 coupled 
NLIEs, exactly one more than for the SU(4) model. 
Their structure should be given by 



4 (£) „oo , 

"EE/ KffH*-v)*B<>\v)% (Al) 

6=lfe=1 ^-oo 

where B^ a \x) = b^\x) + 1 arc the unknown auxiliary 
functions. The free energy is obtained from these func- 
tions via 

4 (a) .OC 1 

/ = -TEE/ V^(y)\nBf\y)^-. (A2) 

a=lj=l- / -° 
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Since the NLIEs must both reproduce the known results 
for the SU(A) model in the limit n — > 1 (fj, — > oo) and 
yield the correct zero-temperature limit, it is possible to 
fix the driving terms and kernel functions. We find 



V {a) {x) = 



Aa 



Ax 2 + a 2 



(A3) 



and the constants 



-2t-n-h/2 
4 1} = c ( P =-2t-n + h/2 



J 2 ) 



-At -2n-h 



J2) _ (2) _ (2) _ (2) 



-At - 2/j, 



J 2 ) 



-At -2n + h 



M _ r (3) _ 

'1 — c 2 — 
,(3) _ J3) _ 



-6t-3/j,-h/2 
-6t-3jj, + h/2 



c[ A) =-8t-Afi. (A4) 
The kernel functions K^ a ^ (x) for a, b — 1,2,3 are similar 



to those of the 5/7(4) model (see Ref. 44 eqs. (33)-(35)), 

>(a,l 



but where the common functions Kj?-,' b \k) are replaced 



by 



= e d-«l*l/» Sinh ^ } - * , 6 . (A5) 
smh(fc/2) 



The remaining kernel functions are 

/oo 
a < 4 )(/fc)e ifea: dfc . (A6) 
-oo 



The set of NLIEs can easily be solved numerically by 
iteration yielding high accuracy over the whole parameter 
range. The validity of the results has been checked by 
comparing our specific heat data to the high-temperature 
expansion rigorously derived in Ref. 45 on the basis of an 
alternative set of NLIEs. Moreover, the results agree in 
the low-temperature limit with CFT and over the whole 
temperature range with numerical TMRG calculations as 
shown in this article. 
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